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1. Introduction 

Quantum Chromodynamics (QCD) phase diagram has been studied extensively, since it is 
related to many fields in physics — the early universe, heavy ion collisions, neutron stars and 
so forth. However, the structure and phase boundary of the QCD phase diagram is under debate 
mainly because of the notorious sign problem at finite chemical potential (/r) in the lattice QCD 
Monte-Carlo simulations. There is a possibility that strong coupling lattice QCD (SC-LQCD) has 
a milder sign problem than the standard lattice QCD simulations, since the effective action can be 
written in terms of hadronic degrees of freedom. In SC-LQCD, the sign problem does not exist 
in the mean field (MF) analysis, and we could sfudy fhe chiral phase fransifion [1-9]. When we 
include flucfuafion effecls, we generally have fhe sign problem as reporfed in monomer-dimer- 
polymer (MDP) simulation [10-12] and auxiliary filed Monfe-Carlo (AFMC) mefhod [13]. The 
sign problem in fhe sfrong coupling limif (SCL) is nol severe in fhese mefhods, and we could 
invesfigafe fhe QCD phase diagram even in fhe finife /r region. 

SCL is fhe opposife limif of fhe confinuum limif, so we need fo include finite coupling effecls lo 
gel insighf info fhe confinuum QCD phase diagram. In fhe MF analyses, we find lhal finite coupling 
[7,8] and Polyakov loop effecls reduce critical lemperalure al /r = 0 [9], and fhe critical poinl (fhe 
endpoinl of fhe firsl order phase fransifion) lemperalure also decreases wilh increasing inverse 
coupling. On fhe olher hand, fhe critical poinl chemical potential evolves in a differenl way in fhe 
cases where we include finite coupling effecls of fhe nexf-fo-leading order (NLO, [7] 

and of fhe nexl-lo-nexl-lo-leading order (NNLO, ^(1 /g^)) [8] in fhe MF Irealmenl. If is necessary 
lo include bolh finite coupling and flucfuafion effecls lo sfudy fhe QCD phase diagram evolulion 
approaching fo fhe confinuum region. Bolh of fhese effecls are evalualed by using reweighing 
mefhod in MDP [12]. Thus if is an imporlanl step lo sfudy fhe QCD phase diagram from direclly 
sampled Monfe-Carlo configurations. In Ihis paper, we conslrucl an effeclive action including finife 
coupling effecls of NLO in AFMC, and evaluate fhe flucfuafion and finite coupling effecls from fhe 
direclly sampled configuralions. 


2. NLO SC-LQCD effective action in Auxiliary field Monte-Carlo method 


In fhe following analysis, we only consider fhe case where fhe spatial lattice spacing a is unify 
(a = 1), and fhe color is 3 (SU(Ac = 3)) in 3-1-1 dimension {d = 3) spacetime. Temporal and spatial 
lattice sizes are indicaled as N-i and L, respeclively. The lattice QCD partition function and aclion 
wilh one species of unrooted slaggered fermion for color SU{Nc) in fhe anisolropic Euclidean 
spacelime are given by 


•^QCD = J ^[X:XiUv]e ^ y^QCD = Sf + Sc 


=5 L L [XxUj,xX,+j - X,+jUlXx\ + ^Lm, , 
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where Xx, Up^ and Up^ represent the quark field, the link variable, and the temporal and spatial 
plaquettes, respectively. The quark chemical potential jJ. is introduced in the form of the temporal 
component of vector potential, rjj x = (— 1 )''°^^ ' is the staggered sign factor, and and 

are mesonic composites. The temporal and spatial physical lattice spacing ratio is introduced as 
/(y) = a/ax. For simplicity, we here assume that the anisotropy parameters in the quark and gluon 
sectors are the same (7 = and the temporal and spatial coupling constants are also the same 
{gx = gs = g)- Under these assumptions, the anisotropy is found to scale the lattice spacing ratio as 
/(y) = y^ and we can define femperafure as T = J^/Nx from fhe mean field argumenfs in SCL [4]. 
We ufilize fhis prescripfion in fhe presenf work. The above laffice QCD acfion has chiral symmefry 
f7(l)v X f7(l)A in fhe chiral limif /tiq —)■ 0 . 

We obfain an NLO effective action affer fhe spatial link infegral as follows [3-8]. 
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“^eff 

^(SCL) 

“^eff 




= + ASW + ^{l/Vd, 1//) , 
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(2.6) 
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( 2 . 8 ) 


where Cx = l/{2N/g'^Y) and Cs = l/{2N/g^Y). We only consider fhe leading order terms of fhe 
large diminufional (l/d) expansion [ 2 ]. 

We now apply fhe exfended Hubaard-Sfrafonovich fransformaion [13] in order fo reduce fhe 
inferacfion terms of quarks info fhe bilinear form. As for fhe spatial NLO ferms, we utilize sequen¬ 
tial bosonizafion. Firsf, fhe bosonized effective acfion for spafial NLO ferms is given as 


„(EHS) _ 
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where = Zk[Yj)^osuk, L/” = 

, and £, = (-1 ^xoH hxj ^hich is related fo 75 in fhe 
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continuum fheory. In Eq. (^^, we divide fhe momenfum region according fo positive ( Ku'^ > 0) 
and negative < 0) modes, and use fhe relation = —Ku\ where u = u + ( 71 , 71 , 71 ) as 
SCL [13]. Second, we use fhe EHS fransformafion fo fhe spafial NEO acfion, sI^loO’ again. 
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We also obfain fhe bosonized acfion terms of temporal NEO terms as 

S™=l3c, £ /(*)[|m*,,|2 + |O*,,|2]+l£[v4(^)Aa7-F0(Ai)Aa4] , (2.11) 

T:.k,f{k)>0 ^ X 
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Aa- = > ( 2 - 12 ) 

j 

Aa+ = Cr'£{^(^x+j + i(Ox-j + i£^)x+j + > (2-13) 

j 

where f{k) = l^jcoskj. (Ox = Zkj{k)>Qe'^''‘^k{T:), = {-^yZk.f(k)>t)e'^''‘^k{T:)- With the SCL 

terms [13], the total effective action after EHS is given as 
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It is clear that the mass and chemical potential are modified by the auxiliary fields. In fhe femporal 
kinetic terms, fhe wave function renormalizafion facfor Zx appears as in fhe MF analysis [7]. 

Nexf, we infegrafe ouf Grassmann variables and femporal link variables analytically, fhen fhe 
effecfive action is written in ferms of hadronic degrees of freedom as 



Me 
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t.U.Kl>0,j 


Iy(i) 12 

\^u \ 
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(2.16) 


where Z(jr) = Yli^xj and ^{x) = Ylifix i ■ We evaluate using fhe recursion formulae [3]. In 
fhe lasf step, we infegrafe ouf auxiliary fields by Monfe-Carlo fechnique af a lime. We generale 
Monle-Carlo configurations and invesligale auxiliary field fluefuafion effeefs in AFMC. 


3. Chiral condensate in NLO with fluctuations 

In order lo undersfand fhe finife coupling effecls from femporal and spafial plaquelles sepa¬ 
rately, we examine fhe phase fransifion al /i = 0 wifh Iwo effective aefions: one is fhe action in 
which we consider temporal NFO and SCF auxiliary fields named as f-NFO, and fhe ofher is fhe 
action in which we fake accounl of spatial NFO and SCF terms named as sp-NFO. 

First, we show the chiral condensate in t-NFO at /r = 0 on a 4^ x 4 lattice in Fig. We here 
regard the chiral radius as the chiral condensate. We find fhal fhe chiral condensate is reduced by 
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chiral cond., AFMC and aniso. MF (t.NLO), (3=0,3 




average phase factor (t.NLO), (3=0, 3 



T 



Figure 1: The chiral condensate (top left), average phase factor (top right), wave function renormalization 
Z (bottom left), and Xn^ (bottom right) as functions of temperature at p = 0 for j3 = 2Nc jg^ = 0, 1, 3 in 
t-NLO. The chiral condensate in AFMC is compared with the MF results, and also with the MF results with 
reduced temperature, 0.93 x T. 


both fluctuation and finite coupling effects. The main finite coupling effect is the mass reduction 
via the wave function renormalization factor Zx as discussed in MF analyses [7]. Since Zx coming 
from the temporal auxiliary fields is larger than unity as shown in Fig. the effective mass nix/Zx 
is suppressed, then the chiral condensate is reduced [7]. By comparison, fluctuations modify 
and play a role to enhance the apparent temperature: the chiral condensate in AFMC agrees well 
with the MF result at enhanced temperature, (pAFMc{T) — /0.93). To our surprise, fluctuation 

effects and finite coupling effects suppress the chiral condensate cumulatively since the both results 
at zero and finite jS in AFMC are well described by 93 percent scaling results in MF. In t-NLO 
simulation, the average phase factor is larger than 0.9 and the wight cancellation is not severe. 
Using the MF knowledge, the repulsive force coming from the auxiliary field, cOt, leads to severe 
sign problem. However, the field (Ot vanishes at zero p and do not affect the sign problem in 
MF. We could expect that the average phase factor is large enough to simulate observables even in 
AFMC at zero p. 

In Fig. 1^ we show the chiral condensate and average phase factor in sp-NLO. The chiral 
condensate takes almost the same value as MF results, while the average phase factor is smaller 
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Figure 2: The chiral condensate and average phase factor as a function of temperature for (5—0,1 with sp. 
NLO. The chiral condensate is almost the same value up to current analysis. Average phase factor is almost 
zero at low temperature, which means the sign problem is so severe. 


than that in t-NLO and is consistent with zero at low temperatures. We need to develop some ways 
to weaken the sign problem in AFMC with spatial NLO terms. 


4. Summary 

We have studied the fluctuation and finite coupling effects on QCD phase diagram in the 
framework of the strong coupling expansion and auxiliary fields Monte-Carlo (AFMC) method. 
We have constructed an effective action of auxiliary fields with both fluctuation and finite coupling 
effects. We apply the Monte-Carlo technique to integrate over auxiliary fields in order to take 
fluctuation effects into account. We have numerically examined the finite coupling effects from 
temporal and spatial plaquettes separately. When we take account of strong coupling limit (SCL) 
and temporal NLO auxiliary fields denoted as t-NLO, both of finite coupling and fluctuation effects 
are found to reduce the chiral condensate. In t-NLO simulations, we find that the average phase 
factor is large enough at zero chemical potential. When we consider SCL and spatial NLO auxiliary 
fields (sp-NLO), the chiral condensate are not modified much in the current analysis. The average 
phase factor is small, and it is necessary to develop a new way to weaken or avoid the sign problem 
in AFMC in order to obtain results in sp-NLO and full NLO simulation. 
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